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Abstract 

Relativistic properties of qq potentail, mass spectra of orbitally and 
radially excited B and D mesons as well as semileptonic decays of B 
mesons to orbitally excited D mesons are discussed in the framework 
of the relativistic quark model based on the quasipotential approach. 

1. Relativistic properties of qq potential. In preceding papers [[I], 
|2|] we have developed the relativistic quark model with the (qq) potential 
consisting of the perturbative one-gluon exchange part and a nonperturbative 
one which is a mixture of the Lorentz scalar and vector confining potentials: 

V(p,q;M) = u a (p)u b (- V )\^a s D^(k)^ al l 

+VZAk)K^ + KonfWWqM-q), (1) 
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where k = p — q, D^ u is the gluon propagator in the Coulomb gauge and T M 
is the effective vector long-range vertex, containing both the Dirac and Pauli 
terms 

r M = 7^ + 7^ a /- r ' ( 2 ) 

ii a> f,(p) are the Dirac bispinors. The parameter (1 + k) can be treated as the 
nonperturbative (long-range) chromomagnetic moment of the quark and n 
as its anomalous part (flavour independent). 

In the nonrelativistic limit the Fourier transform of eq. (HD gives the static 
potential 

V {r) = Vcoui(r) + KonfM + Ko», (3) 



where 



Vcoui(r) = -- — . 

6 r 



= (1 - e)(Ar + B); V* n{ (r) = e(Ar + B), (4) 

here e is the mixing parameter. 

Now assuming that both quarks are heavy enough we evaluate the (v 2 /c 2 ) 
relativistic corrections to the static potential (|3|), (|) paying special attention 
to retardation effects. The Fourier transform of the linear potential At in 
the momentum space looks like: 

A [ d 3 rre~ ik r = -A—^-, k = p - q. (5) 



The natural (though not unique) relativistic extension (dependent only on the 
four-momentum transfer) of expression fl5|) is to substitute (— k 2 ) — > (fcg— k 2 ). 
Now we should choose the procedure of fixing k . On the mass shell due to 
energy conservation we have k = 0. So k may be considered as the measure 
of deviation either from the mass shell or from the energy shell. We choose 
the second possibility and set k equal to e a (p) — e a (q) or to e&(q) — e&(p). 
Then in the symmetrized form we get 



K = ~(e a (p) - e a (q))(e 6 (p) - e 6 (q)), e , 6 (p) = ^Jp 2 + m 2 a b . 

This form is not unique and other possible expressions for /cq are discussed 
in P, M. In favour of choice (Bl) we mention the following arguments. It is 
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well-known || that for the one-photon exchange contribution in QED only 
choice (H) in the Feynman (diagonal) gauge leads to the same correct result 
(the Breit-Fermi Hamiltonian) as the prescription k = in the Coulomb (or 
transverse Landau) gauge. The same is naturally true for the one-gluon ex- 
change contribution in QCD. Moreover as shown in ref. || for any effective 
vector potential generated by a vector exchange and its couplings to con- 
served vector currents (vertices) there is the so-called instantaneous gauge 
which plays the role of the Coulomb gauge. In the instantaneous gauge the 
prescription k = reproduces the same result as the expansion in k 2 . fixed 
by eq. (|6|) in the diagonal gauge used here. The other reason to utilize pre- 
scription (Q) is the reproduction of the correct Dirac limit in this case 0. 

Carring out p 2 /m 2 epansion of the (qq) potential (|I|) in configuration 
space we obtain for the spin-independet part the following expression |], ||: 

V si (r) = V (r) + Vy D (r) + \ ( \ + -\) A[y Coul (r) + (1 + 2*)^ (r)], (6) 

° \ m a m bJ 

where V (r) is given by eqs. (|), (||). The velocity-dependent part VVdM 
can be presented in the form 

VVo(r) = — !— |p 2 H c (r) + i^lv c (r)\ 
m a m h { H ) w 

with 

V bc (r) = -^+{l- £ ) Ar + ( 1 - ^ V ^ = -f 1 - \ Ar > 
V de (r) = --At + Q - I) B; V e (r) = 0. (8) 

Now we are able to test the fulfilment of the exact Barchielli, Brambilla, 
Prosperi (BBP) relations J7|, which follow from the Lorentz invariance of the 
Wilson loop. In our notations these relations look like 

V de -\v bc +\v Q = Q- Ve + ±V e + ~± = (9) 

One can easily find that the functions @ identically satisfy relations @ 
independently of values of the parameters e and k. This is a highly nontrivial 
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result. For the perturbative one-gluon-exchange part of Vyr> our expressions 
for Vb, • • ., V e are the same as in 0, but for the confining (long-range) part 
they are different. The terms with the Laplacian in (|6|) coincide only for 
k = and e = 0, i. e. for purely vector confining interaction without the 
Pauli term in the vertex (^). Our expressions (0) for purely vector (e = 0) 
and purely scalar (e = 1) interactions and for k = coincide with those 
of ref. [U except for the constant B term. Our B term for e — 1 (scalar 
potential) is the same as in JF]. The B term from ref. 0] does not satisfy the 
BBP relations (it gives contribution —B/2 only to Vd e )- Our result for the 
scalar (e = 1) confining potential also differs from the one obtained in ref. [H], 
where the prescription ko = was used and as a result the contribution of 
retardation was lost. The differences between our results and the results 
presented in ref. [§] originate from the use of specific models such as minimal 
area law, flux tube, dual superconductivity and stochastic vacuum. 

The spin- dependent part of our potential [I], §] (for k = — 1) completely 
coincides with the one found in refs. flOl , [7]. The Gromes relation is identi- 



cally fulfilled. Our result supports the conjecture that the long-range con- 
fining forces are dominated by chromoelectric interaction and that the chro- 
momagnetic interaction vanishes. It is also in accord with the flux tube and 
dual superconductivity picture || [10| . 

2. Mass spectra of heavy-light mesons. Many different approaches 
have been used for the calculation of orbital and radial excitations of heavy- 



light mesons [11], [12], However, almost in all of them the expansion in 



inverse powers not only of the heavy quark mass (rriq) but also in inverse 
powers of the light quark mass (m q ) is carried out. The estimates of the light 
quark velocity in these mesons show that the light quark is highly relativistic 
(v/c ~ 0.7 -T- 0.8). Thus the nonrelativistic approximation is not adequate 
for the light quark and one cannot guarantee the numerical accuracy of the 
expansion in inverse powers of the light quark mass. Here we present the 
results of the calculation of the masses of orbitally and radially excited B 
and D mesons without employing the expansion in l/m q (see [jTJJ for details). 
Thus the light quark is treated fully relativistically. Concerning the heavy 
quark we apply the expansion in l/mq up to the first order. Our numerical 
results are presented in Tables 

Let us compare the obtained results with model independent predictions 
of heavy quark effective theory (HQET). Heavy quark symmetry provides 
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Table 1: Mass spectrum of D mesons with the account of 1/mg corrections 
in comparison with other quark model predictions and experimental data. 
All masses are given in GeV. We use the notation (nLj) for meson states, 
where J is the total angular momentum of the meson. 



State 


our 


[11] 


|13| 


experiment ||15|, [16 1 




1.875 


1.88 




1.8645(5) 


is, 


2.009 


2.04 




2.0067(5) 




2.459 


2.50 


2.460 


2.4589(20) 


iPi 


2.414 


2.47 


2.415 


2.4222(18) 


ipi 


2.501 


2.46 


2.585 




IPo 


2.438 


2.40 


2.565 




2S 


2.579 


2.58 






2St 


2.629 


2.64 




2.637(9) ? 



relations between excited states of B and D mesons, such as 

M Bl - M Dl = M Bal - M Dal = M B -M D = M Bs - M Ds = m b - m c , (10) 

where M Bl = (3M Bl +5M B , 2 )/8, M B = (M B + 3M B *)/4 are appropriate spin 
averaged P- and S"-wave states. We get from Tables [l|-|4] the following values 
of mass splittings 

M Bl - M Dl = 3.29 GeV _M Bsl -M Dal = 3.30 GeV 
M B -M D = 3.34 GeV M Bs - M Ds = 3.33 GeV, (11) 



in agreement with (|T0|). There arise also the following relations between 
hyperfme splittings of levels 

TTi Til 

AM B = M B - M B = — -AM D = (M D2 - M Dl ) , (12) 

and the same for B s and D s mesons as well as for P± — P states. Our model 
predictions for these splittings are displayed in Table |^. 

In Tables [l|-^ we compare our relativistic quark model results for heavy- 
light meson masses with the predictions of other quark models of Godfrey 
and Isgur JO]], Isgur [[OJ, Eichten, Hill and Quigg [12[] and experimental data 
[p!5| , [I~6|l . All these quark models use the expansion in inverse powers both of 
the heavy ttlq and light m q quark masses for the Qq interaction potential. In 
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Table 2: Mass spectrum of D s mesons with the account of I/tuq corrections 
in comparison with other quark model predictions and experimental data. 
All masses are given in GeV. 



State 


our 


1 HI 


[12| 


experiment |15|. |16 


15-0 


1.981 


1.98 




1.9685(6) 


ISi 


2.111 


2.13 




2.1124(7) 


IP2 


2.560 


2.59 


2.561 


2.5735(17) 


IPi 


2.515 


2.56 


2.526 


2.535(4) 


IPi 


2.569 


2.55 






lPo 


2.508 


2.48 






2S 


2.670 


2.67 






2St 


2.716 


2.73 







ref. [[TIJ some relativization of the potential has been put in by hand, such 
as relativistic smearing of coordinates and replacing the factors l/m g by 
l/e q (p). However, the resulting potential in this approach accounts only for 
some of the relativistic effects, while the others, which are of the same order 



of magnitude, are missing. The considerations of Refs. [|13|, 12] are closely 
related. The heavy quark expansion is extended to light (u, d, s) quarks and 
the experimental data on P wave masses of K mesons are used to obtain 
predictions for B and D mesons. 

In the paper |13| it is argued that the heavy quark spin P-wave multiplets 
with j = 1/2 (0 + , 1 + ) and j = 3/2 (1 + , 2 + ) in B and D mesons are inverted. 
The 2 + and 1 + states lie about 150 MeV below the 1 + and + states. In 
the limit itlq — > 00, we find the same inversion of these multiplets in our 
model, but the gap between j = 1/2 and j = 3/2 states is smaller (~ 90 
MeV for B and D mesons and ~ 70 MeV for B s and D s mesons) , and 1 / tuq 
corrections reduce this gap further. However, the hyperfine splittings among 
the states in these multiplets turn out to be larger than in [Oj. As a result, 
the states from the multiplets for D, D s and B s mesons overlap in our model, 
however the heavy quark spin averaged centres are still inverted. We obtain 
the following ordering of P states (with masses increasing from left to right): 
B meson — lP^f), 1P 2 , 1P , lP^i); D s meson — 1P , lP^f), 1P 2 , ip(i); 
D and B s mesons - - 1Pi(|), IPq, IP2, 1Pi(|). Thus only for B meson 



we get the purely inverted pattern. Note that the model O] predicts the 
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Table 3: Mass spectrum of B mesons with the account of 1/mg corrections 
in comparison with other quark model predictions and experimental data. 
All masses are given in GeV. 



State 


our 


1 HI 


|13| 




experiment [15, 16] 


15-0 


5.285 


5.31 






5.2792(18) 


ISi 


5.324 


5.37 






5.3248(18) 


m 


5.733 


5.80 


5.715 


5.771 


5.730(9) 


iPi 


5.719 


5.78 


5.700 


5.759 




iPi 


5.757 


5.78 


5.875 






lPo 


5.738 


5.76 


5.870 






2S 


5.883 


5.90 








2St 


5.898 


5.93 






5.90(2) ? 



ordinary ordering of levels. The results of our model agree well with available 
experimental data. 

3. Semileptonic B decays to orbitally excited D mesons. In 
ref. JT7J we have applied the relativistic quark model to the consideration of 
semileptonic decays of B mesons to orbitally excited charmed mesons in the 
leading order of the heavy quark expansion. At the leading order of the heavy 
quark expansion the Isgur-Wise |[8| functions r 3 / 2 and Ty 2 are as follows 



V2 



1 



3 {w + lf/ 2 J (2tt) 3 
d p 



^D(3/2)(P + 



2e fJ 



M D{3/2 )(w + l) 



n/2{w) 



1 



^ 6<2 dp ~^ 
1 



e q + m q 



^b(p) 



(13) 



3y/2(w + iy/ 2 J (2tt)3 
2p 



d 3 P T ( 



2e n 



M, 



D(l/2) 



w + 1) 



-2 — 

q dp 



e + m Q 



^b(p), 



(14) 



where the arrow over 8/ dp indicates that the derivative acts on the wave 
function of the D** meson. The last terms in the square brackets of these 
expressions result from the wave function transformation associated with the 
relativistic rotation of the light quark spin (Wigner rotation) in passing to 
the moving reference frame. These terms are numerically important and 
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Table 4: Mass spectrum of B s mesons with the account of I/tuq corrections 
in comparison with other quark model predictions and experimental data. 
All masses are given in GeV. 



State 


our 


1 HI 


[12| 


experiment |15|. |16 


15-0 


5.375 


5.39 




5.3693(20) 


ISi 


5.412 


5.45 




5.416(4) ? 


m 


5.844 


5.88 


5.861 


5.853(15) ? 


iPi 


5.831 


5.86 


5.849 




iPi 


5.859 


5.86 






IPo 


5.841 


5.83 






2S 


5.971 


6.27 






2St 


5.984 


6.34 







Table 5: Hyperfine splittings of P levels. All values are given in MeV. 



States 


AM D 


— AMfl 

m h ^ 


AM B 


AM Da 


mj, s 


AM Ba 


1P 2 ~ IPi 


45 


14 


14 


45 


14 


13 


1P X - 1P 


63 


20 


19 


61 


19 


18 



lead to the suppression of the T\/2 form factor compared to r^n- Note that if 
we had applied a simplified nonrelativistic quark model |18| these important 



contributions would be missing. Neglecting further the small difference be- 
tween the wave functions ipD(i/2) an d i>D(3/2), the following relation between 
t 3 /2 and ri/2 would have been obtained [|]X 



Tl/ 2 (w) = ^^T 3/2 {w). (15) 

At the point w = 1, where the initial B meson and final D** are at rest, we 
find instead the relation 

tw(i) - n /2 (D ^/^-M^), (16) 

obtained by assuming Vd(3/2) — ipD(i/2) — i>D**- 

In Table ^ we present our numerical results for Tj(l) and its slope in 
comparison with other model predictions ||T9| , |20| , |2l| , p3j [i~Tj |24l . We 
see that most of the above approaches predict close values for the function 
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Table 6: The comparison of our model results for the values of the functions 
Tj at zero recoil of final D** meson and their slopes with other predictions. 



our 


19| 


211 


22| 


m 


Ni, 


HI 


MM 


Ml) 

P3/2 


0.49 


0.41 


0.56 


0.66 




0.54 


0.52 


1.53 




1.5 




2.3 




1.9 




1.5 


1.45 


n /2 (i) 

pi/2 


0.28 


0.41 


0.09 


0.41 


0.35 ±0.08 


0.22 


0.06 


1.04 




1.0 




1.1 




1.4 


2.5 ± 1.0 


0.8 


3 


0.73 



T 3/2(l) and its slope while the results for r 1 / 2 (l) significantly differ from 
each other. This difference is a consequence of a different treatment of the 
relativistic quark dynamics. Nonrelativistic approaches predict t 3 / 2 (1) — 
7-1/2(1) ( see (HI))) while the relativistic treatment leads to t 3 / 2 (1) > 7-1/2(1) 
(see (|i~6D). Our results for the branching ratios of B — > Z?i ]2 (3/2)ez/ decays 
are consistent with available experimental data. 

The work was supported by the Deutsche Forschungsgemeinschaft un- 
der contract Eb 139/1-3 and by the Russian Foundation for Fundamental 
Research under Grant No. 96-02-17171. 
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